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Abstract 

The extended Falicov-Kimball model, with both an on-site hybridization potential and dispersive 
narrow band, is examined within the saddle-point approximation to the Kotliar-Ruckenstein slave 
boson theory. We first set the hybridization potential to zero and find that the phase diagram 
depends strongly upon the orbital structure: for degenerate orbitals, a correlated-insulating state 
is found at sufficiently strong interaction strengths, whereas a finite orbital energy difference can 
lead to discontinuous valence transitions. The obtained phase diagram is very sensitive to the 
presence of a finite hybridization potential. As in Hartree-Fock theory, we find an enhancement of 
the hybridization by the inter-orbital Coulomb repulsion. The more precise treatment of correlation 
effects, however, leads to large deviations from the Hartree-Fock results. In the limit of vanishing 
hybridization an excitonic insulator state is only found when the orbitals are degenerate, which 
restricts this phase to a much smaller parameter space than in other available mean- field theories. 

PACS numbers: 71.30.+h, 71.28.+d, 71.27.+a, 71.10.Fd 
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I. INTRODUCTION 



The Falicov-Kimball model (FKM) was one of the first theoretical attempts to explain 
valence transitions in mixed-valence systems such as SmB 6 and Ce.— In its original form, 
the model describes a spinless fermion system with conduction (c) electrons interacting via 
on-site Coulomb repulsion G with narrow band (/) electrons of orbital energy £/. Within 
a mean-field theory, it was found that varying G or £j can produce both continuous and 
discontinuous changes in the distribution of electrons between these two orbitals, i.e. valence 
transitions. The FKM is nevertheless not a good model of the mixed-valence state as the 
entirely localized nature of the / electrons is unrealistic- The "classical" nature of the / 
electrons was subsequently exploited in the reinterpretation of the FKM as a model of charge 
order in binary alloys. 3 

The central idea behind the FKM, that the inter-orbital Coulomb repulsion G could be 
the origin of discontinuous valence transitions, was revisited by several groups who modified 
the model to account for the quantum nature of the / electrons.-^ The so-called extended 
Falicov-Kimball model (EFKM) allows for the partial derealization of the / electrons due 
to the overlap of the orbital wave functions by including a c-f hybridization potential V. 
Although mean-field impurity models suggested that discontinuous valence transitions were 
possible in the EFKM,-^ weak-coupling studies of the periodic limit found only continuous 
changes in the orbital populations.-^ 

An interesting aspect of the Hartree-Fock (HF) solution of the EFKM with on-site hy- 
bridization is the existence of an excitonic insulator (EI) phase in the limit V — > 0. 5 The 
EI phase is characterized by a spontaneous hybridization between the c and / bands due 
to the presence of a non-zero polarization or excitonic average (c'f). This is interpreted 
as indicating the spontaneous pairing of c electrons with / holes, forming an excitonic con- 
densate. Introduced independently by Keldysh and Kopaev and des Cloizeaux,-^ the EI is 
an important concept in the study of semimetal- semiconductor transitions. A spontaneous 
hybridization between the nested portions of the Fermi surface has proved a particularly 
useful description of the the spin density-wave phase of Cr-based alloys;— a similar scenario 
has also recently been proposed as an explanation for the orbital ordering in LaMnOs.— 
Unambiguous examples of the EI phase remain rare, however, with only two likely candi- 
dates, the alloys Sm 9 La .iS and TmSeo.45Teo.55, thus far discovered.— This indicates that 
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the conditions for an EI phase must be significantly more restrictive than those encountered 
in the usual phenomenological Hamiltonian description. 

The EFKM is the only "standard" tight-binding model of correlated electron systems 
that has been claimed to display an EI phase. 5 It has attracted much attention due to 
the proposal by Portengen et al. that the spontaneous excitonic average in the EFKM 
could be interpreted as evidence of electronic ferroelectricity.— Although a variety of more 
sophisticated treatments'" 5 - or more general mean-field theories^ 1 ^ have failed to find the 
EI phase, the presence of a finite /-electron hopping can stabilize the EI state in the strong- 
coupling regime.— 1 ^ Furthermore, it seems likely that in the EFKM with V 7^ the inter- 
orbital Coulomb interaction will induce a large "excitonic" renormalization of the bare on- 
site hybridization potential.— 1 ^ 

The continuing controversy regarding the EI phase in the EFKM and the larger ques- 
tion about the rarity of EI phases motivates us to study the EFKM using a more advanced 
analytical technique than the weak-coupling methods hitherto employed. A particularly 
useful analytic approach for obtaining the ground state properties of strongly-correlated 
lattice models is the slave boson (SB) mean-field theory developed by Kotliar and Rucken- 
stein. ?0 This is superior to HF mean-field theory as it accounts for the renormalization of 
the quasiparticle weight by the interactions, similar to a Fermi-liquid description. It is also 
of interest to study the SB solution of the EFKM for the possible application to multi-band 
Hubbard models™ in general, the presence of an on-site hybridization makes it impossible 
to apply the usual SB formalism as the atomic Hamiltonian cannot then be written only 
in terms of density operators. Although generalizations of the SB mean-field theory have 
been developed to cope with these difficulties,— the effect of inter-orbital interactions on the 
hybridization is still poorly understood. We can obtain some insight into this situation by 
studying the EFKM with a finite hybridization, as this can be treated within the usual SB 
formulation 

In this paper we examine the EFKM at zero temperature (T = 0) and half-filling using 
the Kotliar- Ruckenstein SB theory. In Sec. [Ill we outline the construction of the mean- field 
SB Hamiltonian as well as reviewing the usual HF solution. In both cases, we consider only 
uniform ground states. The results are presented in Sec. IIHI The solution of the V = 
system (Sec. IIII Aj) is found to be very sensitive to the orbital structure. For degenerate c 
and / orbitals, a Brinkman- Rice-like insulating state is found at sufficiently large interaction 
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strength; for non-degenerate orbitals, discontinuous valence transitions can be found. The 
orbital structure also determines the behaviour of the EFKM with V ^ and tff < 
(Sec. IHI Bp : for degenerate orbitals, the SB solution resembles closely the predictions of the 
HF theory; for non-degenerate orbitals, the more accurate treatment of correlation effects in 
the SB theory produces strong deviations from the HF results. A first-order metal-insulator 
transition is found with increasing interaction strength for tff > ( Sec. IHI C|) . in contrast 
to the second-order transitions found in HF theory. Within the SB treatment of the model, 
the EI phase is only possible when the c and / orbitals are degenerate (Sec. IIIIDI) . We 
conclude in Sec. [TV] with a summary of our results and outlook for further work. 



II. SLAVE-BOSON HAMILTONIAN 



The Hamiltonian for the EFKM is written 

n EFKM = e ^ + E(% e " + e /)^ + vJ2 {4/k + h.c} + g J2 «W (i) 

k k k j 

where c k (Cj) and / k (/■) are the annihilation operators in momentum (real) space for the 
c and / electrons respectively. The c electron dispersion is e k ; the / electron dispersion is 
assumed to be a multiple \tff\ < 1 of the c-electron dispersion. In this work we consider 
hole-like (tff < 0), dispersionless (tff = 0) and electron-like (tff > 0) / electron bands. For 
G = 0, we illustrate the resulting band dispersions as a function of e k in Fig. [TJ For V ^ 0, 
the bands are of mixed c- and /-character: we refer to the upper and lower bands as the 
conduction (C) and valence (V) bands respectively. Note that for tff < the non- interacting 
ground state is insulating for arbitrarily small V, while it is metallic for sufficiently large 
tff > o. 

In the absence of a hybridization potential between the c- and /-electron states we may 
immediately apply the SB technique. In the case when V ^ 0, however, the SB technique 
cannot be straight-forwardly applied as the atomic Hamiltonian (obtained by neglecting all 
kinetic terms in Eq. (TH)) is not diagonal in the electron occupation operators.— To proceed, 
we re-write the system in terms of a diagonal on-site basis 
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FIG. 1: Different band scenarios for the non-interacting EFKM with ej = 0. W is the bandwidth 
of the bare c electron band, (a) c and / electron bands for tff = 0.2, V = 0.0; (b) C and V bands 
for t ff = 0.0, V = 0.05W; (c) C and V bands for t ff = -0.2, V = 0.05W; (d) C and V bands for 
tff = 0.2, V = 0.05W. 



where 



a = 



sgn(e/ + + ) 



V2 



1- 



AV 2 



AV 2 + e 2 



1/2 



AV< 



AV 2 + e 2 



1/2 



(3) 
(4) 



Note that sgn(x + + ) = 1(— 1) for x > 0(x < 0). In terms of the a and b operators we can 
hence re-write the Hamiltonian 

H EFKM = e k {(« 2 + tffP 2 )b[K + (P 2 + t//« 2 )4a k + (1 - t ff )af3 [a[b k + &ta k ] } 

k 

+ J2 ^A a i + E e Ah + G E n » ( 5 ) 

j i 3 

where e b = 6f(3 2 — 2a(3V and e a = tfO? + 2a(3V. 

We adopt the Kotliar-Ruckenstein SB theory by introducing the auxiliary bosonic fields 
£j, s a j, Sbj and dj which respectively destroy the empty, singly-occupied a orbital, singly- 



occupied b orbital and doubly-occupied atomic configurations at site j. The fermionic Hamil- 
tonian is then written in terms of quasi-fermions a and b using the identification 

a j = z aj a j (6) 

bj = Zbjbj (7) 

where 

Za(b)j = (1 - d]dj - si(% S a(%-)" 1/2 (4(a)^i + e ) S a{b) 3 ){^ ~ S l(a)j S b(a)j ~ ^j)^ '• (») 

The physical interpretation of the bosonic fields implies that the following equations are 
satisfied at each site 

1 = + S lj S aj + S bj S bj + d) d 3 ( 9 ) 

= S lj S aj + d ) d 3 ( 10 ) 

W>J = % s 'u ■ (l J (l r ( U ) 

These constraints are respectively enforced by the constraint fields Aj, A aj - and A b j, which 
enter as Lagrangian multipliers. 

A mean-field theory is constructed by replacing the boson and constraint fields by 
spatially-uniform time-invariant fields, etc.. This yields the Hamilto- 

nian 

H SB = ^e k |^ 2 (a 2 + t ff f3 2 )b]j) ]l + z 2 a ((3 2 + t ff a 2 )a\^ v + z a z b (l - t ff )af3 a^k+^k } 

k 

+ ^(e a + A a )4a k + ^(e 5 + A 6 )6^ k 

k k 

+NGd 2 - N\(e 2 + s 2 a + s 2 b + d 2 -l)- NA a (s 2 a + d 2 ) - NA b (s 2 b + d 2 ) (12) 

where 

z a(b) = (l-d 2 - S ^ (b) )- 1/2 (^ 6(o ) + es a(b) )(l - s 2 b(a) - e 2 )- 1 ' 2 (13) 
are the band-renormalization factors. We work throughout at half-filling 

1 = + <*&■>}. ( 14 ) 

as in this limit it may be explicitly demonstrated by extremization of the free energy that 

Z a — %b — Z. 
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The quasi-fermion component of the Hamiltonian Eq. (fl2j) can be straight-forwardly 
diagonalized. For finite c-f hybridization, we have the quasi-fermion C and V bands 

£f (V) = ~ {(1 + t ff )z 2 e k + e a + e b + (-) yj [(1 - t ff )z*e u + e a - ? 6 ] 2 + 4V 2 } (15) 
where 

e a = a 2 (e b + A b ) + (3 2 {e a + A Q ) (16) 

e b = (3 2 {e b + A 6 ) + a 2 {e a + A a ) (17) 

V = a(3(e a + A a -e b -A b ). (18) 

Note the renormalization of the hybridization by the constraint fields: this is the equivalent 
of the excitonic enhancement seen in HF studies. Of particular importance then is the 
so-called excitonic average, defined 

k k 

If A remains finite as V — > 0, the system has an instability towards the EI phase. 

The familiar SB self-consistency conditions are obtained by minimizing the free energy 
with respect to the SB fields while maximizing with respect to the constraint fields, the 
so-called saddle-point approximation. The free energy may be calculated analytically in the 
case of a rectangular density of states (DOS) 

| l/W, U < W/2 
PM= (20) 
I 0, M > W/2 

where the DOS of the bare c and / bands is respectively p c (u>) = p(cu) and pf(u) = 
\t f f\~ l p(t~fj\uj — e/]). Although the details of the self-consistent solutions will change upon 
adopting a more realistic band- structure, the rectangular DOS is convenient for studying 
the generic behaviour of the model. 

A. Hartree-Fock Theory 

For comparison, we briefly discuss the usual HF solution of the EFKM. 6 An effective 
single-particle Hamiltonian may be derived from Eq. ([1]) by decoupling the interaction term 

- Gn B E n / + Gn /E n i- GA E{ c i/i + H - c -} 

3 3 3 3 

-NGn c n f + NGA 2 (21) 



where the HF variational parameters n c , nj and A are respectively the c electron concen- 
tration, the / electron concentration and the excitonic average defined in Eq. f|T9|) . Substi- 
tuting Eq. (1211) into Eq. ([TJ) we obtain the mean-field Hamiltonian 

k k k 

-NGn c n f + NGA 2 (22) 

where <T C = Gnj, e} = e/ + Gn c and = V — GA. The self-consistency equations for 
the HF parameters are easily found by diagonalization of the Hamiltonian Eq. (1221) . see for 
example Ref.j(3. As for the SB results, we calculate the self-consistency equations analytically 
using the rectangular DOS Eq. (I20I) . 

Since the SB parameters are related to the concentration of sites in a given orbital config- 
uration, we can also obtain HF values of these parameters directly from the HF wavefunction 
\^hf)- Explicitly, we have 



d HF = „'52(*HF\n c j n f j \* HF ) 



N 

j 

.2 



= n c (l - n c ) - A (23) 

sIhf = ^ + fnj + ap{c)f 3 + fjc^ HF ) - d 2 HF 

j 

= P 2 + (a 2 - (3 2 )n c - n c (l - n c ) + 2af3A + A 2 (24) 

— Y.^hfIP 2 ^ + « 2 n| - a(3{c]fj + fj Cj }\^ HF ) - d 2 HF 

j 

a 2 + (p 2 - a 2 )n c - n c (l - n c ) - 2apA + A 2 (25) 
1 - s 2 - s 2 - d 2 

1 b aHF b bHF U HF 

n c (l - n c ) - A 2 (26) 



b aHF 



III. RESULTS 



As we work throughout at T = and fixed particle number, we perform the extremization 
upon the ground state energy per site E = (H SB )/N. This is calculated analytically using 
the density of states Eq. f[2"Ul) . The physical values of the SB and constraint fields are then 
obtained by determining the saddle point of the ground state energy, which requires that 
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FIG. 2: (color online) Variation of the slave boson (SB) parameters with G for tff = 0.1, V = 
and various values of e/. The solid lines indicate the ground state solution; the dashed lines give 
the metastable solutions, (a) Band renormalization factor z 2 ; (b) SB field d; (c) SB field s c ; (d) 
constraint field A c — A f . 



we solve the equations 

dE _ dE _dE _dE _dE _ dE _ dE 
de ds a ds b dd d\ dA a dA b 



0. 



(27) 



We employ a multi-dimensional Newton- Raphson technique to solve Eq. fl27|) . To obtain the 
HF results we iterate the self-consistency equations until a desired accuracy is obtained. 



A. V = 0, tff^O 

We begin by examining the EFKM without hybridization. Mathematically, this limit is 
very closely related to the uniform Gutzwiller and SB solutions of the Hubbard model in a 
magnetic field:— >^ identifying the longitudinal magnetic field with e/, the expression for E 
in the two models is of the same form, although the effective bandwidth of the EFKM is 
smaller by a factor (1 + \tff\)/2. This similarity implies that within the SB approximation 
and assuming uniform ground states, the behaviour of the EFKM is a charge analogue of 
the paramagnetic Hubbard model. In particular, there is a localization transition at 6f — 0, 
and at sufficiently small ej^Oa first-order valence transition occurs. This is displayed in 
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our plots of the SB fields in Fig. [21 We do not discuss the HF predictions for the EFKM 
with V = as this only involves the renormalization of the orbital energies, see Eq. f )22|) . 

At tf = there is a second-order transition into a Brinkman- Rice-like correlated-insulator 
(CI) phase at G/W = l + \tff\. This is a localization transition, as the band-renormalization 
factor z 2 vanishes in the CI phase. In the CI phase, every site is singly-occupied with equal 
probability by either a c or / electron, reflected in the limiting values s c = Sf = l/y/2 
in Fig. [21(c) and d = in Fig. 121(b). As implied by the equality of c and / electron populations, 
the difference between the effective c and / energy levels, A c — Af, is zero for all G. 

For any finite e/ the high-G state is an integer valence state with filled c band (/ band) 
for €f > (e/ < 0). This filled-band (FB) state is not localized, so we have z = 1 [Fig. [2(a)]; 
at e/ > it is also characterized by s c = 1 and d = e = Sf = [Fig. E(b,c)]. As G is 
increased, A c — Af decreases to the fixed value e/ — |(1 + in the FB phase: this 

raises the effective / level so that the bottom of the / band just touches the top of the c 
band. 

The transition into the FB state from the metal is an example of a valence transition, 
as it is a transition from a state with mixed valence (i.e. non-zero c and / populations) 
into a state with integer occupation of the orbitals. This transition is qualitatively different 
at small and large |e/|. For small |e/|, the FB state is reached from the low-G metallic 
(M) phase by a first-order transition, which is the charge analogue of the metamagnetic 
transition in the SB treatment of the Hubbard model. Within the region where both FB 
and M solutions of Eq. (T57j) are possible, the stable ground state is defined to be the one with 
the lower energy. Multiple solutions are only found within the region bounded between the 
dashed lines and €f = in Fig. [3j with the metastable state indicated within the brackets 
following the stable state. The FB and M phases have equal energy along the dotted line. As 
|e/| is increased, the boundaries of metastability of the FB and M phases converge together 
at |e/| ~ 0.25W, beyond which we find a second-order transition between the M and FB 
phases. As |e/| is further increased, the line of second order transitions intersects the line 
G = at |e/| = |(1 + t ff )W: for \e f \ > \{l + t ff )W there is no overlap between the c and 
/ bands in the G = limit and so the system is always in the FB state. 



10 



1.4 




' 0.1 0.2 0.3 0.4 0.5 0.6 

|8 f |/W 

FIG. 3: Phase diagram for the EFKM in the G-e f plane. We have = 0.1 and V = 0. At e/ = 0, 
we have a second-order transition between the metal (M) phase and the correlated insulator (CI) 
phase. The limits of metastability of the filled band (FB) and M phases are given by the dotted 
lines; along the dashed line the two phases have equal energy. Where metastable state exists, it 
follows after the stable state in brackets. The solid line for |e/| > 0.25VF gives the second-order 
transition between the M and FB phases. 

B. V^0,t ff <0 

The presence of a finite hybridization potential has a dramatic effect upon the EFKM with 
tff < 0, as the system is then in an insulating state at G = as illustrated in Fig.[]Jb,c). We 
find, however, that the cases £/ = and 6f ^ are distinguished by very different behaviour 
in the V — > limit. 

As shown in Fig. H](a), for tf = and any finite V the band-renormalization factor z 2 
decreases as G is increased, but eventually goes through a minimum before asymptoting to 
unity. For V <C W, z 2 closely follows the band renormalization in the V = system until 
G ~ 0.21F where the curvature changes and a minimum value is subsequently reached at 
G ~ 0.4VF. This has been checked down to V = 10~ 5 W (not shown), which is found to give 
almost identical results to the system with V = 10~ 3 jy. This has an important implication: 
we do not recover the results of Sec. IIII Al in the limit V — > 0, indicating that there is a 
spontaneous hybridization in the V = system. We will discuss this EI phase in more 
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FIG. 4: (color online) Variation of the slave boson (SB) parameters with G for tff = —0.1, ej = 
and various values of V. In (b-d) the solid lines indicate the SB solution while the dotted lines 
indicate the Hartree-Fock solution, (a) Band renormalization factor z 2 ; (b) SB field d; (c) SB field 
s&; (d) change in hybridization Ay. 



detail in Sec. HHP! Although the minimum in z 2 suggests a crossover between two distinct 
regimes, this interpretation has to be used with caution: since the depth of the minimum 
decreases with increasing V, the difference between the low- and high-G regimes becomes 
less pronounced. 

To understand the nature of the low- and high-G regimes we compare in Fig. H^b-d) the 
SB and HF results for d, Sb and the excitonic enhancement of the hybridization AV = V— V. 
For small V, the SB results differ considerably from the HF results for G < 1, although with 
increasing G the results of the two theories converge, which indicates that the physics of the 
high-G regime is HF-like. This is consistent with the G — > oo limit of z 2 , which is unity 
as in a HF theory. This HF-like behaviour is observed for all G at sufficiently large V: for 
V = 0.2W, the SB and HF results for s& and d are almost identical for all G, although there 
is still a noticeable difference in the G-dependence of AV^. This is consistent with the very 
shallow minimum in z 2 observed in Fig. H](a). 

For G < 0.2VF and V = 0.001V4 7 the SB result for d decreases linearly with G, closely 
following the V = results. This is in contrast to the HF theory, which predicts that d 
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differs from its non-interacting value by A ~ exp(— W/G) in this regime. The value of s&, 
however, does not follow the V — results over the same range: this is due to the greater 
sensitivity of sj, to the magnitude of the hybridization, as can be seen by examining the 
evolution of the G = values as V is decreased. Further lowering V we indeed find that the 
V = results are tracked in the low-G regime (not shown). These results, along with the 
variation of z 2 , clearly indicate the importance of correlations beyond HF theory at small 
G and V. 

The different treatment of the on-site Coulomb interaction in the two theories is essen- 
tial to understanding the divergence between the HF and SB results. In the HF theory, 
the Coulomb interaction indirectly affects the concentration of doubly-occupied sites (d 2 ) 
through the renormalization of the band parameters, in this case only the hybridization (at 
€f = the renormalization of the orbital energies is identical). In the SB theory, however, 
the penalty for double occupancy of a site is also explicitly taken into account in the expres- 
sion for the energy by the Gd 2 term. Thus, for G <C W when the Coulomb renormalization 
of the hybridization is exponentially small in both theories, the SB theory nevertheless pre- 
dicts a reduction oc G 2 in d 2 whereas the HF theory has only ~ exp(— 2W/G) reduction. At 
G > W, where the renormalization of the hybridization is large and grows linearly with G 
for both theories, the G-dependence of the SB and HF predictions for d is almost equal. The 
interesting conclusion can then be drawn that within the SB theory the on-site Coulomb re- 
pulsion is less important to the physics than the renormalization of V in the high-G regime. 
That is, the excitonic enhancement of the hybridization gap compensates for the energy 
penalty due to the finite concentration of doubly-occupied sites. 

The G-dependence of the SB fields for ej ^ is qualitatively similar to the €f — 
results for sufficiently large V, but as V is reduced they converge towards the V = results 
of Sec. IIII A[ This can clearly be seen for the V = 0.005W, 6f — 0.05W line in Fig. [5j which 
closely follows the V = curve almost until the valence transition is reached at G ~ 0.7W. 
Although d does not vanish at higher G when V ^ 0, it is heavily suppressed below its value 
in the HF theory. Intermediate between the low- and high-G regimes there is a small range 
of G values for which two solutions exist, implying a discontinuous evolution from the low- 
to the high-G regimes. Unlike the 6f = case [Fig. Hl[c)], for G < 0.7W the G-dependence of 
Sf> is very similar to that of the V = results [dashed line in Fig. Etc)]. The renormalization 
of the hybridization remains small in the low-G regime, but it is much larger and grows 
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FIG. 5: (color online) Variation of the slave boson (SB) parameters with G and V for tff = —0.1 
and ej ^ 0. In all plots the dashed lines indicate unstable solutions; in (b-d) the solid lines indicate 
the SB solution while the dotted lines indicate the Hartree-Fock solution, (a) Band renormalization 
factor z 2 ; (b) SB field d; (c) SB field Sb] (d) change in hybridization AV. 

linearly with G beyond the jump discontinuity. This linear increase is nevertheless slower 
than in the HF theory for all e/ and V. This can be explained by the ^-dependence of the 
hybridization enhancement 

I VI 



AV 



AV 2 + e 2 



=(A«-A 6 ). 



(28) 



In the high-G regime the difference A a — A& grows linearly with the interaction strength. 
Replacing the ^-dependent prefactor by 0.5 we then obtain close correspondence with the HF 
results, as in Fig.H^d). The prefactor in Eq. (1281) . however, arises from the transformation to 
the a-b basis and has important consequences for the SB theory of the EFKM. In particular, 
for finite €f the prefactor vanishes as V — > 0, whereas A a — A& remains finite. This result 
implies that within the SB theory there is no spontaneous hybridization if tf ^ 0, which 
is in strong contradiction to the HF theory where the EI phase is stable at T = for 
|e/| < (1 + tff)W/2^^ Indeed, for Sf ^ 0, the SB and HF theories are only in agreement 
when V ^> |e/|. 

The evolution of the model with decreasing V is shown in the phase diagram Fig. [61 
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FIG. 6: (color online) Phase diagram in the G-V plane for tff = —0.1 and ej = 0.061F. At high G 
we have an excitonically-correlated insulator (ECI), whereas at low G the band-narrowed insulator 
(BNI) is realized. The dashed lines bound the region where both BNI and ECI solutions are found, 
with the metastable state written in brackets. Along the dotted line the BNI and ECI solutions 
have equal energy. The red dot-dash line indicates the location of the minimum in z 2 : within the 
region where both BNI and ECI solutions are found this refers only to the BNI regime. The inset 
shows the phase diagram in the G-tf plane for tff = —0.1 and V = 0.005W. Lines are the same 
as in the main figure. 

We classify the low- and high-G regimes as the band-narrowed insulator (BNI) and the 
excitonically-correlated insulator (ECI) respectively. The former reflects the reduced z 2 < 1 
which is characteristic of the G < W solution when V <C 0.1 W, whereas the latter is 
due to the substantial excitonic renormalization of the hybridization in the high-G regime 
with characteristic V oc G dependence. These designations are most useful where the 
evolution from the BNI to the ECI with increasing G is discontinuous. When there is no 
such discontinuity separating the low- and high-G regimes, the minimum in z 2 can serve as 
an approximate boundary. We emphasize that the BNI and ECI are not distinct phases of 
the EFKM, as both are insulators and there is no order parameter to distinguish between the 
two. Rather, they should be regarded as labeling regions of the phase diagram according 
to the dominant effect of the correlations. We see from the inset that at V — 0.005iy 
multiple solutions are possible within the thin region enclosed by the dashed lines, with 
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the metastable state indicated in brackets. The inset in Fig. [6] can be directly compared 
to Fig. [3j As seen in the main figure, the region of multiple solutions expands to fit the 
V = boundaries as V is decreased, with the BNI (ECI) regimes evolving into the M (FB) 
phase in the V — limit. Note the sensitivity of the region of multiple solutions to a finite 
V: this vanishes completely for V > 0.01PF. We conclude that V ^ strongly suppresses 
any tendency to phase separation in the EFKM. 

C. V^O, t//>0 

The behaviour of the EFKM with V ^ and tff > is qualitatively different to that 
for tff < 0. For sufficiently small V the tff > non-interacting ground state is metallic, 
as shown in Fig. QJd). As G is increased, however, the excitonic renormalization of the 
hybridization eventually opens a gap and the system is then in the ECI phase. In Fig. [7J 
we present representative examples of the different behaviour displayed by the EFKM at 
constant V = 0.05W and e/ = 0. The behaviour of the system for e/ ^ is qualitatively 
identical, although the critical coupling for the metal-insulator transition (MIT) in the SB 
theory is greatly increased at small V due to the much smaller excitonic renormalization of 
the hybridization than at = 0. 

The tff = 0.01 data illustrates the behaviour of the EFKM when the non-interacting 
system is insulating. The system remains in the ECI state for all G, and there are no sig- 
nificant differences between these results and those presented in Sec. IIIIBI More interesting 
is the case tff = 0.1 as there is a second-order MIT at G ~ 0.2W. This is reflected in the 
abrupt change in the first derivative with respect to G of the curves in Fig. [7J Because the 
MIT is driven by the excitonic renormalization of the hybridization, which follows closely 
the HF values in the insulating state, the SB and HF results for the critical coupling agree 
very well. Since the effective hybridization in the SB theory is always smaller than in the 
HF theory, however, the MIT occurs at a slightly higher value of G in the SB theory. 

As tff is further increased the MIT in the SB theory becomes first order. This case is 
represented in Fig. [TJby the tff = 0.5 results. Note that there is a region around G = 0.7W 
where both the ECI and metallic (M) phases are solutions to Eq. ( ]2Tj) . Within the ECI 
phase the SB and HF results are in close agreement for d, Sb and AV^; in the M phase, 
however, the two theories give very different predictions. This is due to the much smaller 
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FIG. 7: (color online) Variation of the slave boson (SB) parameters with G and tff for V = 0.05 
and ej = 0. In all plots the dashed lines indicate unstable solutions; in (b-d) the solid lines indicate 
the SB solution while the dotted lines indicate the Hartree-Fock solution, (a) Band renormalization 
factor z 2 ; (b) SB field d; (c) SB field s&; (d) change in hybridization AV. 

effective hybridization within the SB theory [see Fig. [7(d)] as well as the importance of the 
strong correlations in the M phase. The latter aspect is clearly shown in Fig. [7J(a,b) by the 
close correspondence between z 2 and d in the V = 0.05W and V = systems. Again, this 
is not seen in the G-dependence of s& due to the much greater sensitivity of s& to the value 
of V. Note that as we lower G from the ECI phase and enter into the metastable regime, 
the band gap continuously vanishes as we approach the limit of metastability. Increasing G 
from within the M state, however, the band overlap does not go to zero as the metastable 
limit is approached. 

We present the phase diagram for the EFKM with tff < at V = 0.05W and 6f — 
in Fig. [SI The phase boundaries in both the SB and HF theories are included, the former in 
black and the latter in red; the colour of the labels also refers to the different theories. In 
the HF theory we find only a line of second-order transitions. For the SB theory the MIT 
is of second order for tff < 0.11 and is of first order at higher values of tff, where a region 
with both M and ECI solutions is found. As in the previous phase diagrams the metastable 
phase is given in brackets. The origin of the first-order MIT in the SB theory is related to 
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FIG. 8: (color online) Phase diagram in the G-ttf plane for V = 0.05iy and e/ = 0.0. The red 
line and labels refer to the Hartree-Fock phase diagram, whereas the black lines and labels refer to 
the slave boson theory. The solid line refers to a second-order transitions, the dashed lines to the 
limits of metastability of the excitonically-correlated insulator (ECl) and metal (M) phases. The 
ECl and M phases have equal energy along the dotted line. 

the behaviour of the EFKM in the limit £// = 1. Here the SB solution of the EFKM with 
V 7^ is identical to the paramagnetic SB solution of the Hubbard model in a transverse 
magnetic field, so for sufficiently small V there is a first-order transition into the FB state 
and a region of both M and FB solutions. As t// is reduced from unity, this region narrows 
until the lines of metastability converge at tff ~ 0.11. 

D. Excitonic Insulator 

As noted in Sec. IIIIBI the SB theory only allows an EI phase at e/ = 0, whereas the EI 
phase is a ubiquitous feature of the HF solution. As we see in Fig. [9(a), at e/ = the two 
theories give good agreement for the excitonic average A when G > W and V > 0.1W; for 
G < W and V <C 0.1VF, the SB results predict a smaller |A| as expected from the much 
smaller effective hybridization shown in Fig. 11(d). The predictions of the HF and SB results 
diverge considerably for e/ ^ as displayed in Fig. [9](b). Note that in the HF theory A is 
almost identical for e/ = and e/ = 0.05VF. The very different results of the HF and SB 
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FIG. 9: (color online) (a) Scaling of the excitonic average A with V for the Hartree-Fock (HF) and 
slave boson (SB) theories at tjf = —0.1 and ej = 0. (b) ^-dependence of the excitonic average A 
for the HF and SB theories at t// = -0.1 and G = OAW. 

mean-field approaches suggest separate interpretations of the EI phase. 

In the HF theory the EI phase arises from the formation of an excitonic condensate due 
to the attraction between / holes and c electrons. This excitonic pairing creates an effective 
hybridization between the two bands, as an electron can hybridize from a c orbital into a 
/ orbital via the formation and dissociation of an exciton. At T = the normal (N) state 
is unstable to the EI phase if the effective hybridization is sufficient to open a gap in the 
system. In particular, the EFKM with £// < has an instability towards the EI phase for 
arbitrarily small inter-orbital Coulomb repulsion G. For tff > 0, however, the EI phase is 
only realized when the effective hybridization is large enough to eliminate the band overlap 
and so the N phase is stable up to a finite critical interaction strength. 

To understand the results of the SB theory, we note that the a-b basis introduced in Eq. ([2]) 
rewrites the Hamiltonian in terms of bonding and anti-bonding orbitals at each site. The EI 
phase occurs when these orbitals have mixed c- and /-character in the V — > limit. In the 
non- interacting system this condition is only satisfied for e/ = 0, when the transformation to 
the a-b basis is independent of V, see Eq. (J2J). For 6f ^ 0, in contrast, the transformation to 
the a-b basis is ^/-dependent and the bonding and anti-bonding orbitals continuously evolve 
into the V = atomic orbitals as V — *■ 0. Due to the mixture of the two atomic orbitals 
in the a-b basis, at finite G the Coulomb-induced splitting of the bonding and anti-bonding 
orbitals produces the enhancement of the hybridization given by Eq. ( |28l) . This induced 
hybridization only survives in the V — > limit when e/ = 0, as it is only in this case 
that the atomic orbitals are equally mixed in the a-b basis for all V. As such, we obtain a 
spontaneous hybridization when ej = 0, whereas for ej ^ the system continuously evolves 

19 




FIG. 10: (color online) V = phase diagram in the G-tff plane for ej = 0.0. The red line indicates 
the metal-insulator transition (MIT) between the excitonic insulator (EI) and normal (N) phases 
in Hartree-Fock theory, which are indicated by the red labels. The slave-boson theory predicts a 
first-order MIT between the EI and N phases, with the boundaries of metastability given by the 
dashed lines. The N and EI phases have equal energy along the dotted line. 

into the V = system as clearly evidenced by the phase diagram Fig. 

The 6/ = ground state phase diagram of the EFKM for — l<t//<lis presented 
in Fig. [TUJ Again, we plot the phase diagram for the HF and SB theories in red and black 
respectively. For the HF theory the second-order MIT at tff > found in Sec. IIII Cl remains. 
Although the band gap vanishes continuously as we approach the metallic state from within 
the EI phase, the excitonic average A discontinuously drops to zero at this line and so below 
it the system is in the N state. In contrast to the HF results, the SB theory predicts a 
first-order MIT for all tff > in the limit of vanishing hybridization. As in Sec. IIII CI the 
lower limit of metastability of the EI phase lies just above the MIT in the HF theory, and 
the band gap continuously vanishes as the line of metastability is approached, whereas A 
tends to a finite value. A finite band overlap is always found as we approach the upper limit 
of metastability from within the N phase. Of particular note is the ^-dependence of this 
line: for sufficiently large tff the boundary is given by G = (l+tff)W, which corresponds to 
the limit of metastability of the metallic phase at e/ = + found in Sec. IIII Al At tff ~ 0.1, 
however, there is an abrupt change in this line, with the limit of metastability of the N phase 
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converging toG = 0att// = 0. 



IV. CONCLUSIONS 

In this work we have studied the EFKM within the SB mean-field theory introduced by 
Kotliar and Ruckenstein.— For the system with V 7^ we have compared the predictions 
of SB theory to those of the standard HF approach.- We have found that for the EFKM 
with V = the SB phase diagram displays strong similarities to the SB solution of the 
paramagnetic Hubbard model. In particular, at Sf = we find a transition into a Brinkman- 
Rice-like insulating state as G is increased, whereas at small finite e/ there is a first-order 
valence transition. For finite V a distinction between the EFKM with tff < and tff > 
must be made. In the former case, the model is always in an insulating state and there is 
a considerable renormalization of the hybridization by the Coulomb interaction. For tf = 
the SB and HF results are in good agreement, whereas for e/ 7^ the SB and HF theories 
only coincide when V ^> |e/|. For V <C |e/| the SB theory displays a crossover between a 
state with strong correlations beyond HF-level at G C If and an excitonically-correlated 
state at G > W. At sufficiently large tff > both the HF and SB theories predict a MIT as 
G is increased. In the former this is always of second order; for the latter the MIT is of first 
order for tff greater than some critical value. The presence of first-order transitions in our 
SB treatment suggests that it is worthwhile to include the inter-orbital Coulomb repulsion 
as an important factor in discontinuous valence transitions. 

We have also studied the appearance of an EI phase within the EFKM and our conclusions 
severely constrain the parameter space of the model where such a state is possible. In 
contradiction to the results of HF theory where the EI phase is a ubiquitous feature of the 
T = phase diagram,-^ the SB theory only predicts a spontaneous hybridization when the 
c and / electron atomic orbitals are degenerate, i.e. e/ = 0. When e/ 7^ 0, the effective 
hybridization continuously vanishes as V — > 0. We have explained this difference in terms 
of the importance of the bonding and anti-bonding orbitals in the SB theory. This imposes 
a condition on the atomic structure for the realization of the EI phase. Such a condition is 
absent in the HF theory. 

Our work has only considered uniform ground states of the EFKM. From the rigorous 
solution of the FKM it is expected that the EFKM on a bipartite lattice has an instability 
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towards a density- wave state at e/ = 0. 3 As demonstrated by a number of authors, this 
density-wave state is stable at T = and prevents the EI phase from being realized.— iiL2i 
It is likely that this density-wave phase would also be found within a SB treatment of 
the model. As the condition for the EI phase in the SB theory concerns only the atomic as 
opposed to the band structure, however, the EI phase could still be realized in systems where 
the mean-field density-wave state is unstable, e.g. frustrated lattices. A further extension 
of our work would be to examine the effect of doping away from half-filling. 
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